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Abstract. We update cosmological hot dark matter constraints on neutrinos and 
hadronic axions. Our most restrictive limits use 7-year data from the Wilkinson 
Microwave Anisotropy Probe for the cosmic microwave background anisotropies, the 
halo power spectrum (HPS) from the 7th data release of the Sloan Digital Sky Survey, 
and the Hubble constant from Hubble Space Telescope observations. We find 95% CL 
upper limits of Yl m v < QAA eV (no axions), m a < 0.91 eV (assuming Y^ m f = 0), 
and m v < 0.41 eV and m a < 0.72 eV for two hot dark matter components after 
marginalising over the respective other mass. CMB data alone yield m v < 1-19 eV 
(no axions), while for axions the HPS is crucial for deriving m a constraints. This 
difference can be traced to the fact that for a given hot dark matter fraction axions 
are much more massive than neutrinos. 
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1. Introduction 

Cosmological large-scale structure data allow for precise estimates for the parameters of 
minimal or extended cosmological models. These results have ramifications far beyond 
cosmology itself, notably in the area of neutrino physics. The well-known hot dark 
matter constraints provide neutrino mass limits that directly impact neutrino mass 
searches in single [1] and double [2] beta decay experiments. In a series of papers by our 
collaboration [3-6] and another group [7] this approach was extended to hadronic axions 
where the resulting mass limits are complementary to solar axion searches by the CAST 
experiment [8-10] and the Tokyo axion helioscope [11-13]. In the present work, we use 
the 7-year data release from the Wilkinson Microwave Anisotropy Probe (WMAP) as 
an opportunity to update these results, and also modify along the way several other 
input assumptions as detailed in the main text below. 

Within standard cosmological assumptions, the neutrino plus antineutrino number 
density today, summed over all flavours, is n v ~ 336 cm~ 3 . Currently available 
cosmological data are not yet sensitive enough to resolve the small mass differences 
measured in oscillation experiments, so all neutrinos are treated as having the same 
mass m u , traditionally expressed by the parameter J2 m u = 3m u . For axions, on the 
other hand, the freeze-out temperature and therefore the number density n a depends 
on the axion's interaction rate with pions and nucleons via 

a + 71 <H> 71 + 7T, 

a + N ^ N + 7i, (1.1) 

where the coupling strength is, in turn, proportional to the axion mass m a [4,14,15]. 
Figure 1 shows the relation between m a and n a computed for the thermalisation 
processes (1.1) based on the original calculations of reference [4]. For small m a , the 
number density is also small, so assuming m a = implies n a = which brings us 
back to standard cosmology. Near the hot dark matter limit of m a ~ 1 eV one finds 
a present-day number density of n a ~ 50 cm -3 . Therefore, in the relevant mass range, 
neutrinos are about 6 times more numerous than axions and one expects a hot dark 
matter limit on m a that is roughly twice that on J2 m u, in agreement with what we find 
when we use our full range of input data sets. 

In detail, however, the situation is more subtle. Conventional wisdom says 
that hot dark matter constraints arise primarily from the shape of the measured 
matter power spectrum, since hot dark matter free-streaming suppresses the growth of 
matter perturbations and hence the clustering power on small scales. However, recent 
cosmological data have become so precise that one finds a useful limit on J2 m u of 
order 1 eV already from the cosmic microwave background (CMB) anisotropies alone, 
notably from the increased amplitude of the first acoustic peak in the temperature auto- 
correlation spectrum due to the early integrated Sachs- Wolfe (ISW) effect [16]. The 
same is however not true for axion hot dark matter, since for the same hot dark matter 
fraction, the axion is necessarily some 6 times heavier than the equivalent neutrino. Thus 



Neutrino and axion hot dark matter bounds after WMAP-7 3 
70 i . . . — . . . — 




m a (eV) 

Figure 1. Axion number density function of the axion mass m a assuming the 

hadronic thermalisation processes (1.1) based on the calculations of reference [4]. 

while neutrinos with masses near the hot dark matter limit (m„ < 0.3 eV) essentially 
act like radiation at CMB decoupling and contribute strongly to the early ISW effect, 
the equivalent axion is already nonrelativistic and thus indistinguishable from cold dark 
matter as far as the CMB is concerned. In the latter case, one needs to use the shape 
of the matter power spectrum from smaller-scale data in order to put any sensible 
constraint on m a . 

In order to derive new hot dark matter limits on neutrinos and axions and to explain 
their differences, we begin in section 2 with a description of our cosmological model and 
in section 3 of the data sets used. In section 4 we use standard Bayesian techniques 
to derive credible intervals for m v and m a separately and for a two-component case 
based on different combinations of data sets. We discuss and summarise our findings in 
section 5. 

2. Cosmological model 

We consider a cosmological model with vanishing spatial curvature and adiabatic initial 
conditions, described by eight free parameters, 

= {u cdm ,u h ,H ,T,\n(l0 l0 A s ),n s ,J2 m v, m a}- (2.1) 

Here, Co> c d m = ^cdm^ 2 is the physical cold dark matter density, u h = Q h h 2 the baryon 
density, Hq = h 100 km s _1 Mpc -1 the Hubble parameter, r the optical depth to 
reionisation, A s the amplitude of the primordial scalar power spectrum, and n s its 
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Table 1. Priors and standard values for the cosmological fit parameters considered in 
this work. All priors are uniform (top hat) in the given intervals. 



Parameter 


Standard 


Prior 


^cdm 




0.01-0.99 






0.005-0.1 


h 




0.4-1.0 


T 




0.01-0.8 


ln(10 10 yl s ) 




2.7-4.0 


n s 




0.5-1.5 


X)m„ [eV] 





0-10 


m a [eV] 





0-10 



spectral index. These six parameters represent the simplest parameter set necessary for 
a consistent interpretation of the currently available data. 

In addition, we allow for a nonzero sum of neutrino masses J2 m v and a nonvanishing 
axion mass m a . These extra parameters will be varied one at a time, as well as in 
combination. Their "standard" values are given in table 1, along with the priors for all 
cosmological fit parameters considered here. 

3. Data 

3.1. Cosmic microwave background (CMB) 

We use a compilation of measurements of the CMB temperature and polarisation 
anisotropies from WMAP after seven years of observation [17,18], ACBAR [19], 
BICEP [20], and QuAD [21]. 

3.2. Halo power spectrum (HPS) 

We use the halo power spectrum constructed from the luminous red galaxy (LRG) 
sample of the seventh data release of the Sloan Digital Sky Survey (SDSS-DR7) [22]. 
The full HPS data set consists of 45 data points, covering wavenumbers from /c min = 
0.02 /iMpc" 1 to k max = 0.2 /iMpc -1 . We fit this data set following the procedure 
of reference [22], using a properly smeared power spectrum to model nonlinear mode 
coupling. The smearing procedure requires that we supply a smooth, no-wiggle power 
spectrum, which we construct using the discrete spectral analysis technique introduced 
in reference [23]. 

3.3. Baryon acoustic oscillations (BAO) 

The baryon acoustic oscillation scale has been extracted from SDSS-DR7 [24], which 
provides an angular diameter distance measure at z = 0.275. However, since parameters 
like m v and m a can in principle affect the acoustic scale, care needs to be taken when 
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Table 2. ID marginal 95% upper bounds on m v an d m a for several different choices 
of data sets and models. 



Model 



Fixed m a = 



Fixed m v = 



Data set 



Y. m v [eV] m a [eV] 



CMB only 1.19 

CMB+BAO 0.85 

CMB+HST 0.58 

CMB+HPS 0.61 

CMB+HPS+HST 0.44 



CMB only 

CMB+BAO 

CMB+HST 

CMB+HPS 

CMB+HPS+HST 



Free m v an d m a 



CMB+HPS 0.58 
CMB+HPS+HST 0.41 



No constraint 
No constraint 
No constraint 

1.07 

0.91 



0.82 
0.72 



evaluating the BAO likelihood. We refer the reader to reference [23] for a more detailed 
discussion of this issue. 

3.4- Hubble parameter from the Hubble Space Telescope (HST) 

We adopt the constraint on the Hubble parameter derived from observations with the 
Hubble Space Telescope [25]. 

4. Results 

We use standard Bayesian inference techniques and explore the model parameter space 
with Monte Carlo Markov Chains (MCMC) generated using the publicly available 
CosmoMC package [26]. Our results are summarised in table 2 and figure 2. 

In agreement with several recent papers, the CMB alone provides a fairly robust 
limit on J2 m v With m a held fixed at zero, we find J2 m u < 1-2 eV at 95% C.L. using 
WMAP and other CMB data, while Komatsu et al. find E%< 1.3 eV at 95% C.L. 
from WMAP alone [18]. 

It is noteworthy that for the opposite case when m u is fixed at zero, CMB data 
provide no constraint on m a . At first sight this might seem contradictory, but there is a 
simple explanation. The suppression of small scale power is essentially controlled by the 
fraction of hot to cold dark matter, /hdm = Whdm/w c d m , for both axions and neutrinos. 
This can be seen in figure 4 where the red/solid and green/dash lines have the same 
asymptotic behaviour at large k. However, the number density of axions is six times 
smaller than that of neutrinos. Therefore in order to give the same contribution to Whdm, 
the mass of the axion must be correspondingly larger than the neutrino mass. 

In two examples shown in figures 3 and 4, both the neutrino and the axion hot 




dark matter models have Whdm = 0.013, but the corresponding particle masses are 
J2 m u = 1-2 eV (or m v = 0.4 eV) and m a ~ 2.4 eV respectively. The temperature at 
recombination is approximately 0.2 eV. A neutrino with m„ = 0.4 eV would still be 
semi-relativistic at this time and contribute significantly to the early ISW effect. At 
m a ~ 2.4 eV, however, the axion is nonrelativistic at recombination and behaves more 
like cold dark matter. This disparity between the two hot dark matter candidates can be 
clearly seen in figure 3, where the CMB temperature anisotropy spectrum for the axion 
model is almost indistinguishable from standard ACDM, while its neutrino counterpart 
shows a pronounced early ISW effect. 

The same effect is also manifest in the matter power spectrum shown in figure 4, 
where the epoch of matter-radiation equality — marked by the turning point of the 
spectrum — is different for axions and neutrinos. In the neutrino case, matter-radiation 
equality happens later, leading to a turning point at a much smaller value of k compared 
with its ACDM counterpart. In the axion case, the corresponding matter power 
spectrum traces closely that of the ACDM model up to k ~ 0.02 h Mpc -1 . Beyond this 
point, axion free-streaming sets in and suppresses the power on small scales relative 
ACDM. The amount of suppression in P(k) at large k values is similar in both the 
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Figure 3. CMB temperature anisotropy spectra for various mixtures of cold and hot 
dark matter, where cj c< j m + + w Q = 0.112 is held constant. Bottom: Differences 
to standard ACDM; the shaded/grey area indicates cosmic variance. Light blue/dot- 
dash line: standard ACDM (w cdm = 0.112, u a = 0, ui v = 0). Red/solid line: ACDM+f 
with J2 m v = L2 cV Kdm = 0.099, uj v = 0.013). Green/dash line: ACDM+a with 
m a = 2.4 eV (w c d m = 0.099, u) a = 0.013). Dark blue/dotted line: Extreme axion case 
with m a = 10 eV (cJcdm = 0.0498, u a = 0.0622). 



neutrino and the axion case as expected, since this is governed primarily by /hdm- 

Figures 3 and 4 also show an extreme example of power spectra for m a = 10 eV, 
corresponding to u a = 0.0622, actually exceeding w c d m = 0.0498. Axions with such a 
large mass act essentially as cold dark matter as far as CMB anisotropies are concerned. 

In practice this means that when only CMB data are used, it is difficult to 
distinguish a several-eV axion model from the standard ACDM model, and accordingly 
WMAP provides no useful upper bound on m a . This is true even when additional priors 
are imposed, either in the form of the BAO scale, or the HST measurement of the Hubble 
parameter. Only when low-redshift small-scale data containing shape information such 
as the HPS are added does an upper bound on m a emerge. This is an example of the 
point discussed in reference [23], i.e., that it is crucial to extract the full matter power 
spectrum from a large scale structure survey because it contains important information 
not stored in the geometric BAO scale. 

From our most restrictive data sets and allowing both the neutrino and the axion 
masses to vary simultaneously, we find for the axion mass m a < 0.72 eV (95% C.L.), 
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Figure 4. Same as figure 3, but for the matter power spectrum. 

after marginalising over all other model parameters including Y^ m u- The corresponding 
2D marginalised 68% and 95% contours in the Yli m v~ m a plane are shown in figure 2. 
This new limit on m a is somewhat tighter than our previously published bound of 
m a < 1.02 eV (95% C.L.) using the 5-year WMAP data [6]. Fixing X)m„ to zero, 
the limit on m a becomes less restrictive (m a < 0.91 eV at 95% C.L.), which would be 
applicable if future laboratory experiments were to provide a significant constraint on 
m v . On the other hand, if these experiments were to measure a value of m v close to the 
current cosmological limit, then less room would be left for axion hot dark matter and 
the limit on m a would tighten correspondingly. 

5. Conclusions 

We have provided an updated constraint on axion hot dark matter using new 
cosmological data, most notably CMB data from the WMAP 7-year data release and the 
final SDSS-DR7 LRG data. We have also pointed out a qualitative difference between 
neutrino and axion hot dark matter, in that axion masses in the detectable range are 
large enough that axions are nonrelativistic at recombination. This means they act 
almost like cold dark matter as far as CMB is concerned and that current CMB data 
in itself does not provide a useful limit on the axion mass, even when priors from HST 
or BAO are imposed. In order to properly constrain m a it is necessary to include 
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information on the shape of the matter power spectrum from, for example, the SDSS 
halo power spectrum. 

An interesting question is whether future CMB data will be sensitive to axion hot 
dark matter. The Planck mission has an estimated sensitivity of o~(^2m v ) ~ 0.3-0.5 eV 
[27, 28] . But from figure 3 we can already say that even at high multipoles the effect of 
axion hot dark matter on the CMB anisotropy spectrum is quite small — generally smaller 
than the uncertainty due to cosmic variance. This leads to the inevitable conclusion 
that CMB anisotropy observations will remain poor probes of hot dark matter in the 
several eV range. 

Cosmological bounds on neutrino and axion masses are nicely complementary to 
experimental searches, but cannot replace them. One caveat concerning our axion results 
is that we need to assume that the predicted thermal population was actually produced 
after the QCD epoch. In non-standard cosmologies with low reheating temperature the 
axion population can be severely suppressed and our bounds would not apply [29]. In 
other scenarios a significant cosmic background of low-mass axions can be produced that 
remain relativistic until today and do not form hot dark matter [30]. For neutrinos such 
caveats are less relevant because their thermalisation epoch is well probed by big-bang 
nucleosynthesis and the presence of radiation with roughly the right abundance has been 
confirmed by precision cosmology. 
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